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Ab s tra c t 

An approach for synthesizing buckling results and 
behavior for thin balanced and unbalanced symmetric 
laminates that are subjected to uniform heating or cooling 
and fully restrained against thermal expansion or contrac- 
tion is presented. This approach uses a nondimensional 
analysis for infinitely long, flexurally anisotropic plates 
that are subjected to combined mechanical loads and is 
based on useful nondimensional parameters. In addition, 
stiffness-weighted laminate thermal-expansion parame- 
ters are derived that are used to determine critical temper- 
atures in terms of physically intuitive mechanical 
buckling coefficients, and the effects of membrane orthot- 
ropy and membrane anisotropy are included. Many re- 
sults are presented for some common laminates that are 
intended to facilitate a structural designer’s transition to 
the use of the generic buckling design curves that are pre- 
sented in the paper. Several generic buckling design 
curves are presented that provide physical insight into the 
buckling response in addition to providing useful design 
data. Examples are presented that demonstrate the use of 
the generic design curves. The analysis approach and ge- 
neric results indicate the effects and characteristics of lam- 
inate thermal expansion, membrane orthotropy and 
anisotropy, and flexural orthotropy and anisotropy in a 
very general and unifying manner. 

Introduction 

Buckling behavior of laminated-composite plates 
that are subjected to combined mechanical and thermal 
loads is an important consideration in the preliminary de- 
sign of contemporary, high-performance aerospace vehi- 
cles. The sizing of many structural subcomponents of 
these vehicles is often determined by buckling constraints. 
One subcomponent that is of practical importance in struc- 
tural design is the long rectangular plate. These plates 
commonly appear as elements of stiffened panels that are 
used for wing structures, and as semimonocoque shell 
segments that are used for fuselage structures. Buckling 
results for infinitely long plates are important because 
they often provide a practical estimate of the behavior of 
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finite-length rectangular plates, and they provide infor- 
mation that is useful in explaining the behavior of these 
finite-length plates. Moreover, knowledge of the behav- 
ior of infinitely long plates can provide insight into the 
buckling behavior of more complex structures such as 
stiffened panels. 

An important type of long plate that appears as a 
component of contemporary composite structures is the 
symmetrically laminated plate. In the present paper, the 
term, "symmetrically laminated," refers to plates in 
which every lamina above the plate midplane has a cor- 
responding lamina located at the same distance below the 
plate midplane, with the same thickness, material proper- 
ties, and fiber orientation. Symmetrically laminated 
plates are essentially flat after the manufacturing process 
and exhibit flat prebuckling deformation states, which is 
desirable for many applications. More importantly, the 
amenability of these plates to structural tailoring pro- 
vides symmetrically laminated plates with a significant 
potential for reducing the weight of aerospace vehicles or 
for meeting special performance requirements. Thus, 
understanding the buckling behavior of symmetrically 
laminated plates in a very broad manner is an important 
part of the search for ways to exploit plate orthotropy and 
anisotropy to reduce structural weight or to fulfil a spe- 
cial design requirement. 

For many practical cases, symmetrically laminated 
plates exhibit specially orthotropic behavior. However, 
in some cases, such as thin-walled [±45] s laminates that 
are candidates for spacecraft applications, these plates 
exhibit anisotropy in the form of material-induced cou- 
pling between pure bending and twisting deformations. 
This coupling is referred to herein as flexural anisotropy 
and it generally yields buckling modes that are skewed in 
appearance (see Fig. 1), even when inplane shear loads 
are absent. Symmetrically laminated plates that are un- 
balanced are also being investigated for special-purpose 
uses in aerospace structures. These laminated plates ex- 
hibit anisotropy in the form of material-induced coupling 
between pure inplane dilatation and inplane shear defor- 
mations, in addition to flexural anisotropy. This cou- 
pling is referred to herein as membrane anisotropy and it 
generally yields combined inplane stress states for sim- 
ple loadings like uniform edge compression when in- 
plane displacement constraints are imposed on one or 
more edges of a plate. For example, when the edges of 
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an unbalanced, symmetrically laminated plate, such as a 
[+45/0/90], laminate, are totally restrained against ther- 
mal expansion and contraction, that is caused by uniform 
heating or cooling, inplane shear stresses are developed 
in addition to the usual compressive stresses that are of- 
ten present in balanced laminates. These kinematically 
induced shear stresses may be relatively large in magni- 
tude, compared to the direct compressive stresses, and as 
a result may greatly affect the buckling behavior of the 
plate and yield buckling modes that are skewed in ap- 
pearance. 

The effects of flexural orthotropy and flexural 
anisotropy on the buckling behavior of long rectangular 
plates that are subjected to single and combined mechan- 
ical loading conditions are becoming better understood. 
For example, recent in-depth parametric studies that 
show the effects of flexural orthotropy and flexural 
anisotropy on the buckling behavior of long plates that 
are subjected to compression, shear, pure inplane bend- 
ing, and various combinations of these loads have been 
presented in Refs. 1-3. The results presented in these ref- 
erences indicate that the importance of flexural anisotro- 
py on the buckling resistance of long plates varies with 
the magnitude and type of the combined loading condi- 
tion. Similar results for plates loaded by uniform shear 
and a general linear distribution of axial load across the 
plate width have also been presented in Ref. 4. In a sim- 
ilar manner, the effects of membrane orthotropy and 
membrane anisotropy on the buckling behavior of long 
rectangular plates that are restrained against axial ther- 
mal expansion and contraction and subjected to uniform 
heating or cooling and mechanical loads have been pre- 
sented in Ref. 5. This extensive work has provided a bet- 
ter understanding of the load interaction effects of 
balanced and unbalanced, symmetrically laminated 
plates that are subjected to mechanical loads and re- 
strained against axial thermal expansion and contraction. 

Although an extensive body of work exists that ad- 
dresses the thermal-buckling behavior of plates (see Ref. 
5 for a literature review), a broad understanding of the ef- 
fects of orthotropy and anisotropy on their response has 
not yet been obtained. In particular, the effects of mem- 
brane orthotropy and membrane anisotropy on the buck- 
ling behavior of long rectangular plates that are fully 
restrained against thermal expansion and contraction and 
subjected to uniform heating or cooling are not well un- 
derstood for the large variety of laminated plates that ex- 
ist and the variety of support conditions that are possible. 
One objective of the present paper is to present a more in- 
tuitive buckling analysis for balanced and unbalanced, 
symmetrically laminated plates, that are fully restrained 
against thermal expansion and contraction and subjected 
to uniform heating or cooling. To achieve this goal, the 
buckling analysis is formulated in terms of buckling co- 


efficients for the known, mechanically equivalent loads 
and stiffness-weighted laminate thermal-expansion pa- 
rameters instead of a less intuitive thermal buckling co- 
efficient. Thus, the present study is a continuation or 
extension of the work that is presented in Ref. 5. The 
analysis procedure is based on classical laminated-plate 
theory, which neglects transverse-shear flexibility, and is 
applied to long plates herein. However, the analysis pro- 
cedure is applicable to finite-length plates and to more 
sophisticated plate buckling theories that include effects 
like transverse shear flexibility. 

Two other objectives of the present paper are to 
present a wide range of buckling results in terms of use- 
ful nondimensional design parameters and to provide a 
means for comparing the buckling-response characteris- 
tics of seemingly dissimilar laminated plates. Other ob- 
jectives are to identify the effects of orthotropy and 
anisotropy on the buckling behavior of long symmetri- 
cally laminated plates that are subjected to the same 
loading conditions, and to present some previously un- 
known results. In particular, results are presented for 
plates with the two long edges clamped or simply sup- 
ported and all edges fully restrained against inplane 
movement. Several generic buckling-design curves that 
are applicable to a wide range of laminate constructions 
are presented that use the nondimensional parameters de- 
scribed in Refs. 1-6, along with some other parameters 
that are derived subsequently. Finally, examples are pre- 
sented that demonstrate the use of the generic buckling 
design curves and the analysis procedure. 

Analysis Description 

In preparing generic design charts for buckling of a 
single flat thin plate, a special-purpose analysis is often 
preferred over a general-purpose analysis code, such as a 
finite-element code, because of the cost and effort that is 
usually involved in generating a large number of results 
with a general-purpose code. The results to be presented 
in the proposed paper were obtained by using such a spe- 
cial-purpose buckling analysis that is based on classical 
laminated-plate theory. The analysis details are lengthy; 
hence, only a brief description of the buckling analysis is 
presented herein. First, the buckling analysis for long 
plates that are subjected to a general set of mechanical 
loads is described. Then, the mechanical loads that are 
induced by fully restrained thermal expansion and con- 
traction and that are used in the buckling analysis are de- 
rived, and an expression for the critical temperature 
change is presented that is in terms of the corresponding 
crtitical loading parameter and mechanical-buckling co- 
efficients. 


Symmetrically laminated plates can have many dif- 
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ferent constructions because of the wide variety of mate- 
rial systems, fiber orientations, and stacking sequences 
that can be selected to construct a laminate. A way of 
coping with the large number of choices for laminate 
constructions is to use convenient nondimensional pa- 
rameters in order to understand overall behavioral trends 
and sensitivities of the structural behavior to perturba- 
tions in laminate construction. The buckling analysis 
used in the present paper is based on classical laminated- 
plate theory and the classical Rayleigh-Ritz method, and 
is derived explicitly in terms of the nondimensional pa- 
rameters defined in Refs. 1-6. This approach was moti- 
vated by the need for generic (independent of a specific 
laminate construction) parametric results for composite- 
plate buckling behavior that are expressed in terms of the 
minimum number of independent parameters needed to 
fully characterize the behavior, and that indicate the 
overall trends and sensitivity of the results to changes in 
the parameters. The nondimensional parameters that 
were used to formulate the buckling analysis are given 
by 



(1) 

Q _ Di2 + 2D^ 

(D„D 22 ) l/2 

(2) 

y_ D|6 

(D?,D 22 ) i/4 
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where b is the plate width and X is the half-wave 
length of the buckle pattern of an infinitely long plate 
(see Fig. 1). The subscripted D-terms are the bending 
stiffnesses of classical laminated-plate theory. The 

parameters (X.. and P characterize the flexural orthot- 
ropy, and the parameters y and 5 characterize the flex- 
ural anisotropy. 

The mechanical loading conditions that are includ- 
ed in the buckling analysis are uniform transverse ten- 
sion or compression, uniform shear, and a general linear 
distribution of axial load across the plate width, as de- 
picted in Fig. 1. Typically, an axial stress resultant dis- 
tribution is partitioned into a uniform part and a pure 
bending part. However, this representation is not unique. 
The longitudinal stress resultant N x is partitioned in the 
analysis into a uniform tension or compression part and 
a linearly varying part that corresponds to eccentric in- 
plane bending loads. This partitioning is given by 

= N lc - N b [eo + (e, - eo)t|] (5) 


where N xc denotes the intensity of the constant-valued 
tension or compression part of the load, and the term 
containing N b defines the intensity of the eccentric 

inplane bending load distribution. The symbols ^ and 

e, define the distribution of the inplane bending load, 
and the symbol r| is the nondimensional coordinate 
given by r| = y/b. This particular way of partitioning 
the longitudinal stress resultant was used for conve- 
nience by eliminating the need to calculate the uniform 
and pure bending parts of an axial stress resultant distri- 
bution prior to performing a buckling analysis. 

The analysis is based on a general formulation that 
includes combined destabilizing loads that are propor- 
tional to a positive-valued loading parameter p that is 
increased until buckling occurs, and independent subcrit- 
ical combined loads that remain fixed at a specified load 
level below the value of the buckling load. Herein, the 
term "subcritical load" is defined as any load that does 
not cause buckling to occur. In practice, the subcritical 
loads are applied to a plate prior to, and independent of, 
the destabilizing loads with an intensity below that which 
will cause the plate to buckle. Then, with the subcritical 
loads fixed, the active, destabilizing loads are applied by 
increasing the magnitude of the loading parameter until 
buckling occurs. This approach permits certain types of 
combined-load interaction to be investigated in a direct 
and convenient manner. For example, in analyzing the 
stability of an aircraft fuselage, the nondestabilizing 
transverse tension load in a fuselage panel that is caused 
by cabin pressurization can be considered to remain con- 
stant and, as a result, it can be represented as a passive, 
subcritical load. The combined shear, compression, and 
inplane bending loads that are caused by flight maneu- 
vers can vary and cause buckling and, as a result, they 
can be represented as active, destabilizing loads. 

The distinction between the active, destabilizing 
and passive subcritical loading systems is implemented 
in the buckling analysis by partitioning the prebuckling 
stress resultants as follows 


N, c =-Nj, +N' 2 

(6) 

N y — — N yl + N y2 

(7) 

^xy ^xyl + ^xy2 

(8) 

N b = N bl + N b2 

(9) 


where the stress resultants with the subscript 1 are the 
destabilizing loads, and those with the subscript 2 are 
the subcritical loads. The sign convention used herein 
for positive values of these stress resultants is shown in 
Fig. 1. In particular, positive values of the general linear 
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edge stress distribution parameters N bl , N b2 , and 
e, correspond to compression loads. Negative values 

of N bi and N b2 , or negative values of either Zq or £ h 
yield linearly varying stress distributions that include 
tension. Depictions of a variety of inplane bending load 
distributions are given in Ref. 4. The two normal stress 

resultants of the system of destabilizing loads, N', and 

N yl , are defined to be positive- valued for compression 
loads. This convention results in positive eigenvalues 
being used to indicate instability due to uniform com- 
pression loads. 

The buckling analysis includes several nondimen- 
sional stress resultants associated with Eqs. (6) through 
(9). These dimensionless stress resultants are given by 


N^j b 2 

nxj- K 2 (D n D 22 ) 1 ' 2 

(10) 
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where the subscript j takes on the values of 1 and 2. In 
addition, the destabilizing loads are expressed in terms 

of the loading parameter p in the analysis by 


n c „=L,p 

(14) 

n yl ~ L 2 P 

(15) 

n xyt = p 

(16) 

n bl = L 4 p 

(17) 


where Lj through L 4 are load factors that determine 

the specific form (relative contributions of the load com- 
ponents) of a given system of destabilizing loads. Typi- 
cally, the dominant load factor is assigned a value of 1 
and all others are given as positive or negative fractions. 

Nondimensional buckling coefficients that are used 
herein are given by the values of the dimensionless stress 
resultants of the system of destabilizing loads at the onset 
of buckling; i.e., 
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where the quantities enclosed in the parentheses with 
the subscript “cr” are critical values that correspond to 

buckling and p cr is the magnitude of the loading param- 
eter at buckling. Positive values of the coefficients 
and Ky correspond to uniform compression loads, and 
the coefficient K s corresponds to uniform positive shear. 
The direction of a positive shear stress resultant that acts 
on a plate is shown in Fig. 1. The coefficient K b corre- 
sponds to the specific inplane bending load distribution 

defined by the selected values of the parameters Eq and 
Ex (see Fig. 1). 

The mathematical expression used in the variation- 
al analysis to represent the general off-center and skewed 
buckle pattern is given by 

w N (^,ti) = 1 (A m sinn£ + B m cosn£)<I> ra (r|) (22) 

where £ = xA and r| = y/b are nondimensional coor- 
dinates, w N is the out-of-plane displacement field, and 

A m and B m are the unknown displacement ampli- 
tudes. In accordance with the Rayleigh-Ritz method, 

the basis functions <I> m (r|) are required to satisfy the 
kinematic boundary conditions on the plate edges at r\ = 
0 and 1. For the simply supported plates, the basis 
functions used in the analysis are given by 

<D m (ii) = sin m7CT| (23) 

for values of m = 1, 2, 3, ..., N. Similarly, for the 
clamped plates, the basis functions are given by 

O m (r|) = cos(m-l)7rr| - cos(m+l)rrr| (24) 

For both boundary conditions, the two long edges of a 
plate are free to move in-plane. 

Algebraic equations that govern the buckling be- 
havior of infinitely long plates are obtained by substitut- 
ing the series expansion for the buckling mode given by 
Eq. (22) into a nondimensionalized form of the second 
variation of the total potential energy and then comput- 
ing the integrals appearing in the nondimensional second 
variation in closed form. The resulting equations consti- 
tute a generalized eigenvalue problem that depends on 
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the aspect ratio of the buckle pattern A/b (see Fig. 1) 
and the nondimensional parameters and nondimensiona! 
stress resultants defined herein. The smallest eigenvalue 
of the problem corresponds to buckling and is found by 

specifying a value of X/b and solving the corresponding 
generalized eigenvalue problem for its smallest eigen- 
value. This process is repeated for successive values of 

A/b until the overall smallest eigenvalue is found. 

Results that were obtained from the analysis de- 
scribed herein for uniform compression, uniform shear, 
pure inplane bending (given by £ 0 = -1 and = 1), 
and various combinations of these mechanical loads 
have been compared with other results for isotropic, 
orthotropic, and anisotropic plates that were obtained by 
using other analysis methods. These comparisons are 
discussed in Refs. 1-3, and in every case the results de- 
scribed herein were found to be in good agreement with 
those obtained from other analyses. Likewise, results 
were obtained for isotropic and specially orthotropic 
plates that are subjected to a general linear distribution of 
axial load across the plate width and compared with re- 
sults that were obtained by seven different authors (see 
Ref. 4). In every case, the agreement was good. 


Ereb u c k l in g Stresses.and C ri t i c d Temp e rature.Change 

Uniformly heated or cooled plates that are symmet- 
rically laminated and restrained against thermal expan- 
sion and contraction may develop internal mechanical 
loads that can cause buckling. These induced mechani- 
cal loads enter the analysis through the membrane con- 
stitutive equations. The standard form of these 
membrane constitutive equations for thin plates, that is 
based on classical laminated-plate theory, is found in 
Refs. 7 and 8 and is expressed in terms of membrane 
stiffness coefficients and fictitious thermal stress result- 
ants. In the present study, the membrane constitutive 
equations are used in an inverted form that uses the over- 
all laminate coefficients of thermal expansion and the 
membrane compliance coefficients (see Ref. 8). This 
form of the membrane constitutive equations for sym- 
metrically laminated plates is given by 


v, y 

U>y + v, 3 


a il a i2 a 16 
a i2 a 22 a 26 
a i6 ^6 a 66 


N. 


a. 


N y ) + { a y }0 o 


(25) 


N„ 


a* y 


where u(x,y) and v(x,y) are the prebuckling, inplane dis- 
placements in the x- and y-coordinate directions (see 

Fig. 2), respectively; a x , a y , and a xy are the overall 
laminate coefficients of thermal expansion; the sub- 
scripted a-tcrms are the plate membrane compliance 


coefficients; © 0 is the magnitude of the uniform tem- 
perature change from a predetermined stress- and strain- 
free reference state; and commas followed by a sub- 
script denote partial differentiation with respect to the 
coordinate associated with the subscript. For restrained 
thermal expansion and contraction problems, the plates 
are assumed to be supported and loaded such that the 
prebuckling stress field is uniform, that is, homoge- 
neous. With this assumption, a compatible displacement 
field is obtained directly by integrating Eqs. (25). This 
integration yields 

u(x,y) = |a„ N„ + a, 2 N y + a l6 N ly + a, © 0 jx + g,y + g 2 (26) 

v(x.y) = ^a, 2 N, + a 22 N y + a 26 N xy + a, © 0 jy + g 3 x + g 4 (27) 
and 


g, + g 3 = a 16 N„ + a 26 N y + a 66 N, y + a xy 0 O (28) 

Equations (26)-(28) can be used to determine the 
thermally induced mechanical loadings for several prob- 
lems of practical interest. The problems consist of plates 
restrained against axial thermal expansion or contraction 
(see Ref. 5), plates restrained against transverse thermal 
expansion or contraction (y-coordinate direction), and 
plates restrained against axial and transverse thermal ex- 
pansion and contraction. In the present paper, however, 
only plates that are fully restrained against axial and 
transverse thermal expansion and contraction are consid- 
ered. For this case, all the subcritical loads are zero- val- 
ued and N x (x,y) = -N; is N y (x,y) = -N yl , and N xy (x,y) = 
N xyl (see Fig. 2). All the stress resultants are induced by 
the fully restrained thermal expansion and contraction 
and, when considered together, are destabilizing. The 
displacements given by Eqs. (26) and (27) become 

u(*,y) = (- a„ N;, - a 12 N yl + a I6 N, yl + a, 0 o )x + 

g,y + g 2 (29) 

v(x,y) = (- a, 2 N;, - a* N yl + a M N <yl + a, 0 o )y + 

gjx + g. (30) 

and Eq. (28) becomes 

g . + g 3 = - a i6 N c tl - a* N y , + a* N xyl + a xy © 0 (31) 

Enforcing the restraint condition (displacement bound- 
ary condition) u(0,y) = 0 gives g t = g 2 = 0. Similarly, 
enforcing the restraint condition v(x,0) = 0 gives g 3 = 
g 4 = 0. Enforcing u(a,y) = 0 gives 

" a il “ a i 2 Nyl “*■ a 16 N„ yl + cx, 0 O = 0 (32) 
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( 42 ) 


which yields u(x,y) = 0. Similarly, enforcing v(x,b) = 0 
gives 


^ 12^*1 <*22 N yJ + a 26 N XJf | + ct y 0 O — 0 (33) 

which yields v(x,y) = 0. Substituting g, = g 3 = 0 into Eq. 
(31) gives 

— &i6 N t | — a 26 N y , + N xyl + cc xy 0 O = 0 (34) 

The thermally induced stress resultants for this problem 

are obtained by solving Eqs. (32)-(34) for N x , , N yl , and 
N xyl . The solution is given by 

NI, = (A„a, + A l2 a, + A 1& aj0„ (35) 

N„ = (A l2 a, + A n a y + A^cx „)0 O (36) 

N„ y , = - (A 16 a, + A M a y + A^a J© 0 (37) 

Equations (35)-(37) define a combined loading state that 
is induced by restrained thermal expansion and contrac- 
tion. These equations show that each of the thermally 
induced mechanical loads depends on all three laminate 
coefficients of thermal expansion, and that positive val- 
ues for the stress resultants are possible even for nega- 
tive values of 0 O (uniform cooling) and vice versa. For 
example, a laminate could have a negative value of a x 

and still have a positive value of N x , (axial compres- 
sion). Thus, the signs of 0 O and the parenthetical quan- 
tities in Eqs. (35)-(37) must be considered in 
formulating the buckling problem. 

The buckling problem is posed by first substituting 
Eqs. (35)-(37) into Eqs. (10)-(12), respectively, to obtain 
expressions for the nondimensional stress resultants that 
can be used to characterize the thermally induced me- 
chanical loads. In particular, the nondimensional stress 
resultants are expressed in terms of stiffness-weighted 

laminate thermal-expansion parameters denoted by Otj , 
Ct 2 » and These expressions, with the use of Eqs. 


(14)-(16), are given by 


n I.=^-0.©o = L,p 

(38) 

n yi “ P 

(39) 

n *yi ®3©o — L 3 P 

(40) 

where 


n _ t 2 ( A nCt, + A , 2 a y + A , 6 a xy ) 

1 12VD„D 22 

(41) 


_ t ; (A, 2 a, + A 22 oc y + A a aJ 
2 I2D 22 

t J (A 16 a s + A 26 a, + A w aJ 

1 2(D n D 22 )"‘ (43) 


Equations (10)-(12) and (38)-(40) indicate that N xl , N yI , 

and N xyl are positive- valued when & 2 , and & 3 are 
positive-valued, respectively, and when 0 O is positive- 
valued. Similarly, N xl , N yl , and N xy , are negative- valued 
when , and are positive-valued, respectively, 

and 0 O is negative- valued, or when &i, and & 3 are 
negative-valued, respectively, and 0 O is positive-val- 
ued. Next, Eqs. (38)-(40) are substituted into Eqs. (18)- 
(20) to obtain the relationships between the mechanical- 
buckling coefficients, load factors, and the critical tem- 
perature ®o ; that is, 


K.-(n: i )„ = ^a i 0r = L,p„ (44) 

K,MMc, = ^r« 2 ©o=L 2 p a (45) 

K. s K.L = -^G3©" = L 3 P„ (46) 

where the critical eigenvalue = p a ((3, Y» 8, L„ L 2 , L 3 ) 
for a given set of boundary conditions. For an isotropic 

material, Ct 3 = 0 and <5t, = (5t 2 = a( 1 + v) , where a is the 
coefficient of thermal expansion and v is Poisson’s ratio 
of a homogeneous isotropic material. 

The next step in posing the buckling problem is to 
define the load factors L p L 2 , and L 3 that appear in Eqs. 
(38)-(40) and (44)-(46). It is important to reiterate that 

positive, negative, and zero values for &i© 0 correspond 

to positive, negative, and zero values for N x , , respective- 
ly (see Fig. 2b). Similarly, positive, negative, and zero 

values for & 2 0 o correspond to positive, negative, and 
zero values for N yl , respectively; and positive, negative, 

and zero values for <* 3 ©o correspond to positive, nega- 
tive, and zero values for N xyl , respectively. To define the 

load factors properly, the signs of N xl , N yI , and N xyl must 
be considered. Specifically, the load factors must be de- 
fined such that positive values of L,, L 2 , and L 3 corre- 
spond to positive values of N xl , N yl , and N xyl , 
respectively. Moreover, both positive and negative val- 


6 

American Institute of Aeronautics and Astronautics 


ues of © () must be considered. These requirements lead 
to six cases that must be considered in formulating the 

buckling analysis; that is, the cases for which &i© 0 > 0 , 

a,0 o <o, a, =0 with a 2 0 o >o, a,=o with a 2 © 0 <o, 

a s = (X 2 = 0 with a 3 0 o >O, and <2, = <5 2 = 0 with 

(Xj0 o < 0 . The buckling analysis for each of these cases 
is presented subsequently. 

For the case when 3j© 0 > 0 , N*, > 0 and = 1 is 
appropriate (axial compression). The values for the 
other two load factors that are needed to completely 
define the prebuckling stress state are obtained by divid- 
ing Eqs. (19) and (20) by Eq. (18), with L, = 1, or by 
dividing Eqs. (39) and (40) by Eq. (38). This step yields 

L = ty* ( Pu l _ + A^PCy + A^cc^ f D n ) (Al\ 

2 K[d 22 J a,~A u cc t + A t3 <x y + A l6 a v [D J 1 J 

L = Hqrt f Pul _ ^ _ A lfe cc, + A^cty + A 66 a, y ( p l \ 

3 K\D n J a x A ll tx, + A ll a y + A l6 a x ^D a J ( ) 

For an isotropic material, these expressions reduce to L 2 
= 1 and L 3 = 0. With L, = 1 and L 2 and L 3 defined by 
Eqs. (47) and (48), the relationship between the critical 

value of the mechanical loading parameter p CT and the 

critical temperature is determined by Eq. (44); that 
is, 

K .=^ a . 0 r=po ( 49 ) 

where = p CT (p, y, 6, L ls L 2 , L 3 ) for a given set of 
boundary conditions. It is important to point out that 

Eq. (49) yields positive and negative values for ©^ for 

positive and negative values of respectively. More- 
over, it is important to reiterate that the relationship 

between p CT and the corresponding mechanical buckling 
coefficients K^, Ky, and K s is given by Eqs. (49), (45), 
and (46), respectively. 

For the case when 6i0 o < 0 , N*, < 0 and L, = -1 
is appropriate (axial tension). Like for the previous case, 
the values for the other two load factors that are needed 
to completely define the prebuckling stress state are ob- 
tained by dividing Eqs. (19) and (20) by Eq. (18), but 
with L t = -1, or by dividing Eqs. (39) and (40) by Eq. 
(38). This step yields 

L N >.fD.r_ Q >_ A'i a . + A n a > + D,, \ l/2 

n :.\dJ Q, A „a, + A„a, + A „a,,\ D„ j 


L N «>'fP..V' 4 = gj_ A,.a. + A M a,^-A^a„ (D, l y ,4 
’ n;,\dJ a, A 11 a, + A,,a, + A 1 .a„tDj (5I) 

With Lj = -1 and L 2 and L 3 defined by Eqs. (50) and 
(51), the relationship between the critical value of the 

mechanical loading parameter p cr and the critical tem- 
perature ©o is again determined by Eq. (44); that is, 

K . = ^fi. 0 " = -Pa (52a) 

In contrast to the previous case, Eq. (52a) yields positive 

values for ©o for negative values of and vice versa. 
For laminates with L 3 = 0 (balanced laminates) and 

L 2 ^ 0 (transverse tension), no destabilizing compres- 
sion or shear loads are present and buckling cannot 
occur because the plate is in a state of biaxial tension. In 
contrast, when L 3 = 0 and L 2 > 0, a plate is subjected to 
axial tension and transverse compression loads. Figures 
29-31 of Ref. 1 indicate that an infinitely long plate 
buckles as a wide column for this type of loading, and 
that the buckling coefficient Ky = 1 and 4 for simply 
supported and clamped plates, respectively. With Ky 
known, Eqs. (45) and (52a). give 


K _12b 2 ft0 c_ K y 

l0 °- LT 


(52b) 


For a general symmetric laminate, the possibility 
exists that Gj = 0, which implies that N‘, = 0. For this 


case, Lj = 0 is appropriate and the sign of 6t 2 0 o must be 
considered in defining the nonzero load factors. In par- 
ticular, for & 2 ©o > 0 , N yl > 0 and L, = 0 and L 2 = 1 are 
appropriate (transverse compression). The value for the 
load factor L 3 that is needed to completely define the pre- 
buckling stress state is obtained by dividing Eq. (20) by 
Eq. (19), with L 2 = 1, or by dividing Eqs. (40) by Eq. 
(39). This step yields 


'dJ 

w _ ft, _ A ,.a, + A ls a,+ A^oc 

r D„\ 


1 ft, A,, a, + A u a, + A,,a„ 

Id..; 


(53) 


The relationship between L 3 and the mechanical-buck- 
ling coefficients Ky and K s is shown in Figs. 24-27 of 
Ref. 1. With L t = 0, L 2 =1, and L 3 defined by Eq. (53), 
the relationship between the critical value of the 

mechanical loading parameter Per and the critical tem- 
perature ©r is determined by Eq. (45); that is, 
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(54a) 


(46); that is, 


K > = -pr a 3©o = Pa 


(57) 


It is important to point out that Eq. (54a) yields positive 
and negative values for ©o for positive and negative 

values of G 2 , respectively. For laminates with L 3 = 0 
(balanced laminates) a plate is subjected to only trans- 
verse compression. Thus, an infinitely long plate buck- 
les as a wide column for this type of loading, and the 
buckling coefficient Ky = 1 and 4 for simply supported 
and clamped plates, respectively. With Ky known, Eq. 
(45) gives 

~^x®2©o = K, (54b) 

For the case when (5, = 0 and & 2 ©o < 0 , N yl < 0 
and L, = 0 and L 2 = -1 are appropriate (transverse ten- 
sion). Like for the previous case, the value for the load 
factor L 3 that is needed to completely define the prebuck- 
ling stress state is obtained by dividing Eq. (20) by Eq. 
(19), with L 2 = -1, or by dividing Eqs. (40) by Eq. (39). 
This step yields 

L — _ ^ ,yl f jjii l — ^ _ A lt (X, + 

rj - “ ^ = A .A + A 22 o, + A »<x*A j (55) 

Like the previous case, the relationship between L 3 and 
the mechanical-buckling coefficients Ky and K s is shown 
in Figs. 24-27 of Ref. 1. With L, = 0, L 2 = -1, and L 3 
defined by Eq. (55), the relationship between the critical 

value of the mechanical loading parameter p a and the 

critical temperature ©IT is determined by Eq. (45); that 
is, 

K, = ~^pr®2©o = - P<j (56) 

For this case, Eq. (56) yields positive values for ©o for 

negative values of G 2 t and vice versa. For laminates 
with L 3 = 0 (balanced laminates), no destabilizing com- 
pression or shear loads are present and buckling cannot 
occur because the plate is in a state of uniaxial tension. 

For the case with Ctj = <5t 2 = 0 and G 3 © 0 > 0 , 

N‘, = N yl = 0 and N syl > 0 1 which implies that L, = L 2 = 
0, and that L 3 = 1 is appropriate (positive shear loading). 
With L, = 0, L 2 = 0, and L 3 = 1 , the relationship between 
the critical value of the mechanical loading parameter 

p a and the critical temperature ©" is determined by Eq. 


Again, it is important to point out that Eq. (57) yields 
positive and negative values for ©o for positive and neg- 
ative values of G 3 , respectively. In addition, values of 
the mechanical-buckling coefficient K s for several lami- 
nates are given in Ref. 1. 

The final case to consider is when ft, = ft 2 = 0 and 

ft 3 © 0 <0. For this case, N‘,=N yl =0 and N xy , <0, 
which implies that Lj = L 2 = 0, and that L 3 = - 1 is appro- 
priate (negative shear loading). With = 0, L 2 = 0, and 
L 3 = -1, the relationship between the critical value of the 

mechanical loading parameter and the critical tem- 
perature ©o is again determined by Eq. (46); that is, 

K. = ^jy<S 3 ©o = - Per (58) 

In contrast to the previous case, Eq. (58) yields positive 

values for ©" for negative values of & 3 > and vice versa. 

It is important to mention that the approach used 
herein to define the prebuckling stress state and the crit- 
ical temperature ©^ also applies for a more sophisticated 
plate theory, like a first-order transverse-shear deforma- 
tion theory. For this theory, p„ would depend also upon 
additional nondimensional parameters that characterize 
the transverse-shear flexibility, Thus, the only differ- 
ence in the results for the two plate bending theories is 

the actual value of p a that is used in Eqs. (44)-(46), for a 

given problem. It is also important to point out that p CT 
for an infinitely long plate does not depend on the param- 
eter ct_ . This fact has been shown in Refs. 1-4. 

Results for Common Laminates and Discussion 

Results are presented in this section and in Figs. 3- 
13, that illustrate behavioral trends for several common 
symmetrically laminated plates that are fully restrained 
against thermal expansion and contraction and subjected 
to uniform heating or cooling. In particular, results are 
presented for several common balanced, symmetric lam- 
inates that are made of IM7/5260 graphite-bismaleimide 
material (see Table 1); that is, [(±45/0/90) m ] s quasi-iso- 
tropic laminates, [(±45/0 2 ) m ] s and [(±45/90 2 )J s quasi- 
orthotropic laminates, and [(±0) m ] s angle-ply laminates, 
where a positive value of the lamina fiber angle 0 is 
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shown in Fig. 6. The [(±45/0 2 ) m ] t and [(±45/90 2 )J s lam- 
inates are described as quasi-orthotropic because of the 
presence of a relatively small amount of flexural anisot- 
ropy. Results are also presented for [±0/O/9O] s laminates 
with angle plies and for a quasi-isotropic laminate whose 
principal material-coordinate frame is rotated by an an- 
gle 0; that is, [(±45/0/90) + 0] s . In addition, results are 
presented for similar unbalanced, symmetric laminates; 
primarily, [(+45 2 /0/90)J s , [(+45 2 /0 2 )J s , [(+45 2 /90 2 )J s , 
[(+0) 2m ] s > and [+0 2 /O/9 O] s laminates that exhibit a signif- 
icant degree of membrane anisotropy in addition to flex- 
ural anisotropy. All of the results are based on classical 
laminated-plate theory and the nominal ply thickness 
used in the calculations was 0.005 in. 

Results are presented in Figs. 3-4 that show the 
stiffness-weighted laminate thermal expansion parame- 
ter and load factor ratio as a function of the 

number of laminate plies for several balanced and unbal- 
anced, symmetric laminates. In particular, the black sol- 
id lines in the figures are for the [(±45/0/90)J s , 
[(±45/0 2 ) m ] s , and [(±45/90 2 )J s balanced, symmetric lam- 
inates and three of the gray dashed lines are for the 
[(+45/0/90U, [(+45 2 /0 2 )J s , and [(+45/90,) J s unbal- 
anced, symmetric laminates. The difference in these two 
groups of laminates is the presence of membrane anisot- 
ropy that is caused by orienting all the 45-degree plies in 
the same direction. Results are also presented for 
[(±45/10 2 ) m ] s and [(±45/30 2 )J s unbalanced, symmetric 
laminates with two additional gray dashed lines in Figs. 
3-4. Similarly, results are presented in Fig. 5 that show 

the shear-load-factor ratio as a function of the 

number of laminate plies for the same unbalanced, sym- 
metric laminates. The solid black lines in Fig, 5 are for 
the [(+45 2 /0/90)J s> [(+45/0 2 ) m ] s , and [(+45 2 /90 2 )J s lam- 
inates and the dashed gray lines are for the [(±45/10 2 ) m ] s 
and [(±45/30 2 ) m ] s laminates. The reason for presenting 

( and <VCt, in Figs* 4 and 5, respectively, is be- 
cause these quantities determine the mechanical-buck- 
ling load factors L 2 (see Eqs. (47) and (50)) and L 3 (see 
Eqs. (48) and (51)), respectively, that are used to deter- 
mine the critical temperature change when <5t ( is nonze- 
ro. 

The results in Fig. 3 indicate that the values of Gt, 
are all positive. Thus, these laminates are loaded by axial 
compression when uniformly heated and by axial tension 
when uniformly cooled (see Eqs. (38) and (10)). The re- 
sults also show equal values of G, for the [(±45/0 2 )J s 
and (+45 2 /0 2 ) m ] s laminates, for the [(±45/0/90)J s and 
[(+45/0/90)^ laminates, and for the [(±45/90 2 )J 3 and 
[(+45 2 /90 2 )J s laminates, respectively. The largest val- 


ues of Gi are exhibited by the [(±45/90 2 ) m ] s and 
[(+45 2 /90 2 )J s laminates and the smallest by the 
[(±45/0 2 )J # and [(+45 2 /0 2 )J s laminates. The results also 

show, for the most part, relatively small variations in <2| 
with the number of laminate plies, with the largest varia- 
tions exhibited by the [(±45/0/90) J s and [(+45 2 /0/90) m ] s 
laminates. 

The results in Fig. 4 indicate that the values of 
are all positive, which means that the values for 

& 2 are positive because the values of are all positive. 
Thus, the restrained thermal expansion and contraction 
induces a uniform compressive stress in the y-coordinate 
direction for uniform heating and a uniform tension 
stress for uniform cooling (see Eqs. (39) and (11)). The 

results also show equal values of Q.Jd ] for the 
[(±45/0 2 ) m ] s and [(+45 2 /0 2 ) m ] s laminates, for the 
[(+45/0/90) m ] s and [(+45 2 /0/90) m ] s laminates, and for 
the [(±45/90 2 ) m ] s and [(+45 2 /90 2 ) m ] s laminates, re- 
spectively. The largest values of G/ft, are exhibited by 
the [(±45/0 2 ) m ] s and [(+45 2 /0 2 ) m ] s laminates (<V&i = 
3.2) and the smallest by the [(+45/90 2 ) m ] s and 

[(+45 2 /90 2 ) m ] s laminates (<V<3 m = 0.28). The results 
also show, to a large extent, relatively small variations in 

<V<X, with the number of laminate plies. However, the 
results in Fig. 4 also show a substantial, monotonically 

increasing variation in G.JQ. l with the number of lami- 
nate plies for the [(±45/0 2 ) ra ] s , [(+45 2 /0 2 ) m ] s , and 
[(+45/10 2 ) m ] s laminates. 

The results in Fig. 5 also indicate that the values of 
are all positive, which means that the values for 

are positive because the values of <*i are all positive. 
Thus, the restrained thermal expansion and contraction 
induces a uniform positive shear stress for uniform heat- 
ing and a uniform negative shear stress for uniform cool- 
ing (see Eqs. (40) and (12)). The results also show nearly 

equal values of <V<5t, for the [(±45/10 2 ) m ] s and 
[(+45 2 /90 2 ) m ] s laminates for m > 2. The largest values 
of &y/(X l are exhibited by the [(+45 2 /0 2 ) m ] s laminates 
(a/G, = 0.53) and the smallest by the [(±45/1 0 2 ) m ] s 
laminates (<V<*i = 0.15). The results in Fig. 5 also show 
relatively small variations in with the number of 
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laminate plies. The largest variation in G/G, with the 
number of laminate plies is exhibited by the 

[(+45 2 /0 2 ) m ] s laminates. 

Three curves are presented in Fig. 6-7 that show the 
stiffness-weighted laminate thermal-expansion parame- 
ter G, and load factor ratio G/G, as a function of the fi- 
ber angle, 0° < 0 < 90°, for [(±0) m ] s , [±0/O/9O] s , 

[(±45/0/90) + 0] s balanced laminates and for [(+0) 2m ] s 
and [+0^0/90], unbalanced laminates. Specifically, the 

black solid lines in the figures are for the [(±0) m ] s and 

[(+6) 2 Js laminates and are independent of the value of 
the stacking sequence number, m. The dashed lines are 
for the [±0/0/90], and [+02/0/90], laminates, and the sol- 
id gray curves are for the [(±45/0/90) + 0], laminates. 
Likewise, results are presented in Fig. 8 that show the 

shear-load-factor ratio G/G, as a function of the fiber 

angle 0 for the [(+0)2m] s and [+0 2 /O/9O] s unbalanced 
laminates. In Fig. 8, the solid and dashed black curves 

are for the [(+6)2™] s and [+0 2 /O/9O] s laminates, respec- 
tively. 

The results in Fig. 6 indicate that all of the lami- 
nates exhibit G, > 0, for all values of 0 considered. 
Thus, these laminates are also loaded by axial compres- 
sion when uniformly heated and by axial tension when 

uniformly cooled. The largest and smallest values of G, , 
and greatest variations, are exhibited by the [(±0) m ] s 
and [(+®)2m]s laminates. Both the [(±0) m ] s and 
[(+0)2mL and the [±0/O/9O] s and [+02/0/90], laminates 
exhibit a monotonic increase in G, with increasing val- 
ues of 0. The values of G, for the [(±45/0/90) + 0] s lam- 
inates decrease monotonically with increasing values of 
0 up to 45° and then increase monotonically. Moreover, 
the solid gray curve for these laminates is symmetic 
about the vertical line 0 = 45°. 

The results in Fig. 7 indicate that the values of 
G/G, are all positive, which means that the values for 

G 2 are positive because the values of G, are all positive. 
Thus, the restrained thermal expansion and contraction 
induces a uniform transverse compressive stress N yl for 
uniform heating and a uniform transverse tension stress 
for uniform cooling. The largest and smallest values of 

G/G,, and greatest variations, are exhibited by the 


[(±9) ml, and [(+0) 2m] S laminates. Specifically, the 

largest and smallest values are given by G/G, = 14.39 

and 0.07, respectively. Both the [(±0)m] s and 

[(+0)2m] s and the [±0/O/9O] s and [+0,/O/9O] s laminates 

exhibit a monotonic decrease in G/G, with increasing 
values of 0, approaching a value of zero. The variation 

in the values of G/G, with fiber angle for the 
[(±45/0/90) + 0] s laminates is benign compared to the 
other laminates. 

The results in Fig. 8 also indicate that the values of 
G/G i are all positive, which means that the values for 

G 3 are positive because the values of G, are all positive. 
Thus, the restrained thermal expansion and contraction 
induces a uniform positive shear stress for uniform heat- 
ing and a uniform negative shear stress for uniform cool- 
ing. Overall, the results show that the values of G/G, 

are greater for the [(+0)2m] s laminates than for the 
[+0 2 /O/9O] s laminates, with respect to the fiber angle 0. 

In addition, the results show larger variations in G/G, 
for both laminates for values of the fiber angle less than 

approximately 45°. The largest value of G/G,, exhibited 

by the [(+0)2mL laminates, is given approximately by 

G/G, = 1. 16 at approximately 0 = 21°. The largest value 

of G/G, exhibited by the [+0 2 /O/9O] s laminates is given 

approximately by G/G, = 0.43 at approximately 0 = 26°, 
For all of the unbalanced laminates considered in 
Figs. 3-8, biaxial compression and positive shear loading 
exists when the laminates are subjected to uniform heat- 
ing. For uniform cooling, biaxial tension and negative 
shear loading exists in the laminates. Moreover, a biax- 
ial compression stress state exists in the balanced lami- 
nates subjected to uniform heating, and a biaxial tension 
stress state exists in the balanced laminates subjected to 
uniform cooling, which means that buckling cannot oc- 
cur. Furthermore, it was determined that elastic buckling 
does not occur for the unbalanced laminates considered 
in Figs. 3-8 that are subjected to uniform cooling. 

Results are presented in Figs. 9 and 10 that show 

12b 2 

the critical temperature change "^T®o for uniform 
heating as a function of the number of laminate plies for 
the [(±45/0/90)J s , [(±45/0 2 )J s , and [(±45/90 2 ) m ] s 
balanced laminates and the [(+45 2 /0/90) m ] s , 
[(+45 2 /0 2 ) m ] t , [(+45 2 /90 2 ) m ] s , [(±45/1 0 2 ) m ] s , and 
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[(±45/30 2 ) m ] s unbalanced laminates. Similarly, results 
are presented in Figs. 11 and 12 that show the critical 
temperature change for uniform heating as a function of 

the Fiber angle 0 for the [±9] s , [( + ®)2m] s> and [(±0) m h 
(m > 5) laminates, and for the [±0/O/9O] s , [+0 2 /O/9O] s , 
and [(±45/0/90) + 0] s laminates, respectively. The re- 
sults in Figs. 9 and 10 correspond to results for plates 
with simply supported and clamped edges, respectively. 

The results were obtained by computing the ratios <xj(x x 


and<VG, First. Next, the load factors L 2 and L 3 were de- 
termined by using Eqs. (47) and (48), respectively. Then, 

the critical value of the loading parameter p CT was deter- 
mined for each of the laminates and the critical tempera- 
ture change was obtained by using Eq. (49). Equations 

(47)-(49) were used because &t©o > 0 for these lami- 
nates when subjected to uniform heating. For these lam- 
inates, the critical value of the loading parameter 
depends on the flexural boundary conditions, the plate 
flexural orthotropy and flexural anisotropy, and the val- 
ues of the load factors L 2 and L 3 , which depend on the 
stiffness-weighted laminate thermal-expansion parame- 
ters. The stiffness-weighted laminate thermal-expansion 
parameters depend on the plate membrane orthotropy, 
membrane anisotropy, and the stiffnesses associated 
with pure bending action. 


The results in Figs. 9 and 10 indicate that the criti- 
1 2b 2 

cal temperature is highly dependent on the ar- 

rangement of the 45-deg, 0-deg, and 90-deg plies. 
Generally, as the number of plies increases, the magni- 
tude of the critical temperature change decreases for 
about half of the laminates and increases for the other 
half. Overall, the plates with clamped edges are more 
buckling resistant than the corresponding plates with 
simply supported edges, as expected. The clamped and 


simply supported [(±45/90 2 ) m ] s laminates require the 
most heating to cause buckling and the simply supported 


[(±45/0 2 ) m ] s , [(+45 2 /0 2 ) m ] s , and [(±45/10 2 ) m ] s lami- 


nates and the clamped [(+45 2 /0 2 ) m ] s laminates require 
the least amount of heating. 

The results in Fig. 1 1 indicate that the critical tem- 
perature change -^yp-©o for the [±0] s , [(+9)2mL, and 

[(-®) m ] s (m > 5) laminates is highly dependent on the 
Fiber angle 0. In particular, the results for the clamped 

and simply supported [±0] s and [(±6) m ] s (m > 5) lami- 
nates show substantial increase in the critical tempera- 


ture change with increasing 0 for values up to 
approximately 55° and 62°, respectively, followed by a 

signiFicant decrease. In addition, the [(±0) m ] s (m > 5) 
laminates require more heating to cause buckling than 
the corresponding [±0] s laminates. The results for the 

[(+6)2m] s laminates show a monotonic increase in the 
critical temperature change with 0, but the temperature 
changes are always less than or equal to those for the cor- 
responding [±0] s and [(±0) m ] s (m > 5) laminates. Thus, 

the [(+9)2m] s laminates generally require less heating to 
cause buckling than the other corresponding laminates 
shown in Fig. 11. Overall, the plated with clamped edges 
are much more buckling resistant than the corresponding 
plates with simply supported edges, as expected. 

The results in Fig. 12 indicate that the critical tem- 
perature change ^r©^ for the [±0/O/9O] s , [+02/0/9013, 

and [(±45/0/90) + 0] s laminates is also highly dependent 
on the Fiber angle 0. SpeciFically, the results for the 
[±0/O/9O] s laminates exhibit a trend similar to the results 
for the balanced laminates shown in Fig. 11; that is, the 
clamped and simply supported [±0/O/9O] s laminates ex- 
hibit a substantial increase in the critical temperature 
change with increasing 0 for values up to approximately 
57° and 67°, respectively, followed by a signiFicant de- 
crease. Like the unbalanced laminates of Fig. 1 1, the re- 
sults for the [+02/O/9O] s laminates show a monotonic 
increase in the critical temperature change with 0, but the 
temperature changes are always less than or equal to 
those for the corresponding [±0/O/9O] s laminates. Thus, 
the [+0 2 /O/9O] s laminates require, for the most part, less 
heating to cause buckling than the corresponding 
[±0/O/9O] s laminates. Generally, the [(±45/0/90) + 0] s 
laminates require more uniform heating to cause buck- 
ling than the other corresponding laminates shown in 
Fig. 12, and exhibit greater values of the critical temper- 
ature change over a larger range of 0. Moreover, the crit- 
ical temperature change for the [(±45/0/90) + 0] s 
laminates exhibits the smallest variations with Fiber an- 
gle for the laminates shown in Fig. 12. Like for the lam- 
inates of Fig. 1 1, the laminates in Fig. 12 with clamped 
edges are much more buckling resistant than the corre- 
sponding plates with simply supported edges, as expect- 
ed. 

The data presented in Figs. 9- 1 2 are given in a com- 
pact form and are applicable to an inFinite range of plate 
width-to-thickness ratios, b/t. Once the number of lami- 
nate plies is selected, the critical temperature change ©q 
can be found as a function of the plate width b for a giv- 
en laminate family. However, it is important to keep in 
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mind the limitations of classical laminated-plate bending 
theory as the plate width-to-thickness ratio b/t becomes 
smaller than a value of approximately 20. Results of this 
type are useful in structural design and are shown for the 

[(±45/0/90) m ] s quasi-isotropic laminates in Fig. 13. 
Two sets of curves are shown for various values of the 
stacking sequence number m in Fig. 13. Moreover, the 
results are shown for a maximum temperature-change 
magnitude of 300 °F. This maximum temperature- 
change magnitude was selected to ensure a regime of 
heating in which the material behavior is linear. The sol- 
id lines in these figures are for plates with simply sup- 
ported edges and the dashed lines are for plates with 
clamped edges. As indicated by Eq. (49), the magnitude 

of shown in the figure decreases proportionally to the 
inverse square of the plate width b. 

Generic Results and Examples 

The simplicity of the equations presented herein 
that are like Eq. (49) in form suggested that generic 
buckling design charts, similar to those presented in 
Refs. 1-5, could be obtained for the thermal buckling 
problem of the present study. Specifically, generic de- 
sign data can be obtained from charts that show the stiff- 
ness-weighted laminate thermal-expansion parameters 

t <X 2 , and G 3 ; the flexural orthotropy parameter p; and 
the flexural anisotropy parameters y and 5 as a function 
of material system and laminate stacking sequence, and 
from charts that show the critical value of the loading pa- 
rameter p CT , or buckling interaction curves, as a function 
of flexural boundary conditions, p, y, 5, L,, L 2 and L 3 . 
Results of this type illustrate the key aspects of the be- 
havior and show overall trends and sensitivity of the be- 
havior to changes in the parameters. Several figures that 
illustrate the utility of this type of design data are pre- 
sented subsequently. 

Values of the Nondimensional Parameters 

Values of the nondimensional parameters that are 
used herein are presented in this section for 

[(±45/0/90) m ] s quasi-isotropic, [(±0) m ] s angle-ply, 

and [(+9)2mL unidirectional off-axis laminates. Nine 
different contemporary material systems are used. These 
material systems include boron-aluminum, S-glass-ep- 
oxy, atypical boron-epoxy, AS4/3501-6 graphite-epoxy, 
AS4/3502 graphite-epoxy, IM7/5260 graphite-bismale- 
imide, Kevlar 49-epoxy, IM7/PETI-5, and P-100/3502 
pitch-epoxy materials. The mechanical properties of 
these material systems are presented in Table 1 and the 


nominal ply thickness is 0.005 in. 

Parameters for [(±45/0/90) m ] s laminates. Values 

of G, and for [(±45/0/90) m ] s quasi-isotropic 

laminates are presented in Figs. 14 and 15, respectively, 
for the nine material systems. The results show a wide 

variation in <5t t with material system for the laminate 

considered. Moreover, is positive for all the material 
systems except the P-100/3502 pitch-epoxy material. 

Similarly, G/G, is positive for all the material systems 
and 1 < G/G, < 1.2. These values for <V&, indicate 

that & 2 is also positive for all the material systems ex- 
cept the P-100/3502 pitch-epoxy material. Thus, the 
laminates made from all the material systems except the 
P-100/3502 pitch-epoxy material are loaded by biaxial 
compression when uniformly heated and by biaxial ten- 
sion when uniformly cooled, and as a result can buckle 
only for uniform heating. The laminates made from the 
P-100/3502 pitch-epoxy material are loaded by biaxial 
tension when uniformly heated and by biaxial compres- 
sion when uniformly cooled. These laminates must be 
cooled to buckle. The results also show that some mate- 
rial systems yield a slightly larger variation in G, and 

<V&, with the number of laminate plies. Overall, Figs. 
14 and 15 show that there is a wide range of possibilities 
available for tailoring the thermal buckling characteris- 
tics of a laminate family. 

The remaining parameters needed for a buckling 

analysis of the long [(±45/0/90) m ] s quasi-isotropic lam- 
inates are the the nondimensional flexural orthotropy pa- 
rameter P and the flexural anisotropy parameters y and 
8. Values for these parameters are not presented herein. 
The parameters are given in Ref. 5 as a function of the 
stacking-sequence number m for the 9 material systems, 
along with a discussion of their characteristics. The re- 
sults in Ref. 5 for these parameters show a series of 
curves that approach P = 1 , y = 0, and 5 = 0 from above 
as the number of plies increases. A homogeneous, iso- 
tropic material has values of P = 1 , y = 0, and 6 = 0. 
Thus, the results for P, y, and 8 in Ref. 5 give, to some 
extent, a quantitative measure of quasi-isotropy. 

Parameters for [(±^)mh and laminates. 

Values of and & 3 for [(±0) m ] s balanced, angle- 

ply laminates and [(+0) 2m ] s unbalanced, unidirectional 
off-axis laminates composed of the 9 material systems 
are presented in Figs. 16-18, respectively. The values for 

Ct, and & 2 arc identical for both laminates and indepen- 
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dent of the stacking sequence number m. The results 

show a very wide variation in 6| , , and & 3 with ma- * 

terial system and with fiber angle 0. The largest varia- 
tions are exhibited by , followed by . Moreover, 

and & 2 are positive for all the material systems except 
the P-100/3502 pitch-epoxy and Kevlar 49-epoxy mate- 
rials, which are negative for several values of the fiber 

angle 0. Similarly, & 3 is positive for all the material sys- 
tems except the S-g lass-epoxy and boron-epoxy materi- 
als, which are negative for 0° < 0 < 90°. Thus, depending 
on the material system and fiber angle, a laminate may be 
loaded by various combinations of tension, compression, 
and shear and may buckle when subjected to uniform 
heating, cooling, or both. Like the results for the quasi- 
isotropic laminates, the results in Figs. 16-18 show that 
there is a very wide range of possibilities available for 
tailoring the thermal buckling characteristics of a lami- 
nate family. 

Values of the nondimensional orthotropy parame- 
ter P and the nondimensional anisotropy parameters y 
and 6 for the [(±0) m ] s laminates made of the same 9 ma- 
terial systems have also been presented in Ref. 5. The 
values of P in Ref. 5 for the [(±0) m ] s laminates are iden- 
tical to the corresponding results for the [(+0) 2m ] s unbal- 
anced, unidirectional off-axis laminates and are 
independent of the stacking-sequence number m. The 
results in Ref. 5 show a series of curves for P that vary 
dramatically with the fiber angle 0 and the material sys- 
tem. The largest values of, and greatest variations in, P 
are generally exhibited by the laminates made of the P- 
100/3502 pitch-epoxy material. In contrast, the smallest 
values of, and least variations in, P are generally exhib- 
ited by the laminates made of the boron-aluminum mate- 
rial. The values of y and 5 for the f(±0) m ] s laminates are 
strongly dependent on the stacking-sequence number m. 
The results in Ref. 5 for the flexural anisotropy parame- 
ters y and 5 are for m = 1, which corresponds to the 
highest degree of flexural anisotropy for these laminates. 
The results in Ref. 5 show a large effect of the fiber angle 
0 and the material system on the degree of flexural 
anisotropy for the [±0] s laminates. Like the parameter p, 
the largest values of, and greatest variations in, y and 5 
are generally exhibited by the laminates made of the P- 
100/3502 pitch-epoxy material, and the smallest values 
of, and least variations in, y and 6 are generally exhib- 
ited by the laminates made of the boron-aluminum mate- 
rial. 

Values of the nondimensional anisotropy parame- 
ters y and 5 for the [(+0) 2 J S laminates made of the same 
9 material systems are presented in Figs. 19 and 20, re- 
spectively. The results for these laminates are indepen- 


dent of the stacking sequence number m and exhibit 
practically the same trends as the results for correspond- 
ing [±0] s laminates that are given in Ref. 5. The values 
of y and 5 for the [(+0) 2m ] s laminates, however, are gen- 
erally larger than the corresponding values for the 
[(±0) m ] s laminates. This fact is illustrated in Fig. 21 
which shows a comparison of results for [(±0) 6 ] s , [±0] 5 , 
and [(+0) 2 J S laminates made of IM7/5260 graphite-bis- 
maleimide material. The results in this figure indicate 
that the [(+0) 2m ] s unidirectional off-axis and [±0] s lami- 
nates exhibit similar variations with 0, but the [(+0) 2m ] s 
laminates exhibit a substantially higher degree of flexur- 
al anisotropy. In addition, the results for the [(±0) 6 ] s lam- 
inates indicate that the flexural anisotropy is negligible 
compared to that of the other laminates. 

Buckling Coefficients and Critical Temperature Change 

The values of the nondimensional parameters p, y 
and 5 given in Ref. 5 and presented in the present paper, 
or from similar figures for other laminates, can be used 
to determine the buckling coefficients for plates subject- 
ed to uniform, combined mechanical loads. More specif- 
ically, the critical loading parameter P« and the buckling 
coefficient Ky, and K s depend on the parameters P, y, 
and 5 for classical laminated-plate theory (see Ref. 1) 
and the load factors L } , L 2 , and L 3 that define the relative 
proportions of the inplane loads. For plates subjected to 
restained thermal expansion and contraction and uniform 
heating or cooling, the relative proportions of the in- 
plane loads that are induced by the restrained deforma- 
tion are defined by the stiffness-weighted laminate 

thermal-expansion parameters 0t 2 , and & 3 (e.g., see 
Eqs. (47) and (48)). Thus, generic buckling design charts 
for plates subjected to mechanical loads can be used di- 
rectly to obtain critical temperature changes for the ther- 
mal buckling problem considered herein, and can be 
applied to a vast range of laminate constructions that in- 
clude hybrid laminates made of several materials. Sev- 
eral of these generic buckling design charts that are ap- 
plicable to balanced and unbalanced, symmetric 
laminates are presented subsequently. 

Charts for balanced laminates. Balanced laminates 
that are subjected to uniform heating or cooling and fully 
restrained against thermal expansion and contraction de- 
velop, at most, a biaxial stress state (G 3 = 0). For the case 

when G, and G 2 are both positive valued, a laminate is 
loaded by uniform biaxial compression when heated uni- 
formly. In contrast, for the case when G, > 0 and & 2 < 
0, a laminate is loaded by uniform axial compression and 
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transverse tension (N yl < 0) when heated uniformly. A 

similar stress state exist when when <X t < 0 and <X 2 > 0 
and the laminate is uniformly cooled. However, for all 
three of these loading conditions; - 1, L 3 = 0, and L 2 
is given by Eq. (47). The critical temperature change 
1 2b z 

~^ 2 ®o is obtained by dividing both sides of Eq. (49) by 

ft, , once p CT (|3, y, 5, L„ L 2 , L 3 ) is known. The appropri- 
ate value for is obtained from generic buckling de- 
sign charts by noting that Eqs. (44) and (45) give K, =p„ 

and Ky = L 2 p CT , respectively, and that L 2 = Ky/K^. Thus, 
the critical value of the loading parameter for balanced 
laminates can be obtained from traditional Kj-Ky buck- 
ling interaction curves that are given as a function of the 
nondimensional parameters P, y, 8, and the boundary 
conditions. 

For the case when and are both positive val- 
ued and the laminate is subjected to uniform cooling, or 

the case when and & 2 are both negative valued and 
the laminate is subjected to uniform heating, a state of bi- 
axial tension exists in the laminate. For this loading con- 
dition, elastic buckling is not possible. For the case when 

> 0 and &2 < 0, a laminate is loaded by uniform axial 
tension and transverse compression when cooled uni- 
formly. A similar stress state exist when when Gt, < 0 

and & 2 > 0 and the laminate is uniformly heated. For 
both of these cases; Lj = - 1 , L 3 = 0, and L 2 is given by Eq. 
(50). More importantly, because the only destabilizing 
load is a transverse compression load, a wide-column 
buckling mode is the only possibility. Thus, the critical 
1 2b 2 

temperature change is obtained by dividing 

both sides of Eq. (52b) by G, . For this special case, Ky 
= 1 and Ky = 4 for simply supported and clamped plates, 
respectively, and K* = - Ky/L 2 . 

For the case when <5fc, = 0 and 6 2 < 0, a laminate is 
loaded by uniform transverse tension when heated uni- 
formly. A similar stress state exist when when Ct, = 0 

and > 0 and the laminate is uniformly cooled. For 
this loading condition, elastic buckling is also not possi- 
ble. For the cases where = 0 and <X 2 < 0 and the lam- 
inate is cooled uniformly, and where Qj = 0 and > 0 
and the laminate is heated uniformly, a state of transverse 
compression exists. A wide-column buckling mode is 
the only possibility. The critical temperature change 


1 2b 

is obtained by dividing both sides of Eq. (54b) 

by . For this special case, Ky = 1 and K y = 4 for sim- 
ply supported and clamped plates, respectively. 

Examples of traditional K x -Ky buckling interaction 

curves that can be used to obtain p CT (P, y, 8, L„ L : , L 3 ) 
described in the previous paragraphs are presented in 
Fig. 22 for specially orthotropic plates, that have negligi- 
ble flexural anisotropy (y = 8 = 0). In particular, two sets 
of generic buckling curves are shown in Fig. 22 for 
plates that are subjected to combined axial compression 
and transverse tension or compression loads. Results are 
shown in the figure for six values of the orthotropy pa- 
rameter P that cover a very wide range of laminate con- 
structions. The solid and dashed curves correspond to 
results for clamped and simply supported plates, respec- 
tively, and the curves for P = 1 correspond to results for 
isotropic plates. Similar curves that show the effects of 
flexural anisotropy are presented in Ref. 1, along with a 
discussion of the behavioral characteristics. For the K^- 
Ky buckling interaction curves shown in Fig. 22, the load 
factor L 2 appears implicitly as the slope of a line emanat- 
ing from the origin of the graph. The buckling coeffi- 
cients that correspond to a given value of L 2 are 
determined by the point of intersection of the line and the 
appropriate buckling interaction curve. These points of 
intersection are depicted in Fig. 22 for the simply sup- 
ported and clamped plates by the open and filled circular 
symbols, respectively. Points of intersection that occur 
on the horizontal portion of the curves, given by Ky = 1 
and 4 for the simply supported and clamped plates, re- 
spectively, correspond to wide-column buckling modes. 
Similarly, results for negative values of K^, although not 
shown in Fig. 22, correspond horizontal extrapolations 
of the lines Ky = 1 and 4, and also correspond to wide- 
column buckling modes. With the value of L 2 deter- 
mined from Eq. (47), unique values of K^ and Ky are de- 
termined from a figure like Fig. 22, where = K*. 

Results that are directly analogous to those present- 
ed in Fig. 22 are presented in Figs. 23 and 24 in which 

the buckling coefficient K x , or equivalently p CT , is given 
explicitly in terms of the load factor L 2 . An interesting 
characteristic of Fig. 23 is that the family of curves for 
the simply supported plates (dashed) and the clamped 
plates (solid) terminate at points on the gray, solid curves 

given by K X = J- and K x = ^, respectively. The 

curves K, = J- and K x = y- correspond to the horizon- 
l 2 

tal portion of the curves in Fig. 22 that are given by Ky = 
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I and K y = 4 for the simply supported plates (dashed) and 
the clamped plates (solid), respectively, and represent 
the wide-column buckling modes. Like in Fig. 22, neg- 
ative values of K x (not shown) correspond to wide-col- 
umn buckling modes. 

Typically, when studying the behavior of plates 
that possesses flexural anisotropy, one likes to know the 
importance of the anisotropy. Generic curves for simply 
supported plates that can be used in conjunction with 
Figs. 23 and 24 to obtain this information are presented 
in Figs. 25 and 26. The curves indicate the effects of 
plate flexural anisotropy on the ratio of the buckling co- 
efficients that include and neglect this anisotropy, re- 
spectively, as a function of (3 and L 2 . Buckling 
coefficients that correspond to the neglection of flexural 

anisotropy are denoted by KJ T . 5 . 0 . Six groups of 
curves are shown in each figure that correspond to equal 
values of the flexural anisotropy parameters. Moreover, 
three curves are contained in each group that correspond 
to different values of L 2 . 


Charts for unbalanced laminates. Unbalanced lam- 
inates that are subjected to uniform heating or cooling 
and fully restrained against thermal expansion and con- 
traction can generally develop a prebuckling stress state 
that consists of combinations of axial tension or com- 
pression, transverse tension or compression, and positive 
or negative shear. The specific form of the prebuckling 
stress state is determined by the values of the parameters 

G i , 0 - 2 , and ; which determine the values of the load 
factors L p L 2 , and L 3 . The specific values of L p L 2> and 
L 3 are obtained by following the same logic presented 
previously for the balanced laminates and discussed dur- 
ing the presentation of Eqs. (47)-(58). Once L p L 2 , and 

L 3 are determined, pcr(P> Y» 5, L„ L 2 , L 3 ) can be deter- 


mined from generic results that give either or K^-Ky- 
K s buckling interaction surfaces (or level curves) as a 
function of L p L 2> L 3 , p, y, 5, and the boundary condi- 

1 2 

tions. Then, the critical temperature change is 

determined by using the appropriate one of Eqs. (49), 
(52a), (54a), (56), (57), and (58). 


Examples of generic buckling curves are shown in 
Fig. 27 for plates with negligible flexural anisotropy (y = 
6 = 0) and subjected to combined axial compression, 
transverse tension or compression, and shear. These re- 
sults are applicable to unbalanced, symmetric laminates 
with negligible flexural anisotropy that develop biaxial 
and shear loads when restrained from thermal expansion 


and contraction and can be used to obtain p a described 


in the previous paragraph. The results shown in Fig. 27 
represent the traditional buckling interaction curves and 
are for a value of P = l, which corresponds to results for 
isotropic plates. The load factors L 2 and L 3 are included 
in the results in an implicit manner. The solid and dashed 
curves correspond to results for clamped and simply sup- 
ported plates, respectively, and are level curves of the 
corresponding K x -K y -K s buckling interaction surface. 
Each curve in the figure corresponds to a different mag- 
nitude of the shear loading. Moreover, positive and neg- 
ative values of the shear-buckling coefficient K s 
correspond to positive (see Fig. 2b) and negative direc- 
tions of the shear loading. The buckling coefficients that 
correspond to given values of the load factors L 2 and L 3 
are determined by the point of intersection of the line em- 
anating from the origin of the graph and the appropriate 
buckling interaction curve. These points of intersection 
are depicted in Fig. 27 for the simply supported and 
clamped plates by the open and filled circular symbols, 
respectively. Points of intersection that occur on the hor- 
izontal portion of the curves, given by Ky = 1 and 4 for 
the simply supported and clamped plates, respectively, 
correspond to wide-column buckling modes. With the 
values ofL„ L 2 , and L 3 determined from the values of the 

parameters G, , <X 2 , and % , unique values of K x , K y , and 

K s are determined from a figure like Fig. 27, and p a is 
determined from Eqs. (44)-(46). It is worth mentioning 
that very few results of this type are available in the tech- 
nical literature. 

Results that are directly analogous to those present- 
ed in Fig. 27 are presented in Figs. 28-33 for simply sup- 
ported plates in which the buckling coefficient K^, or 

equivalently p tr , is given explicitly in terms of the load 
factors L 2 and L 3 , with Lj = 1 . Thus, the results in these 
figures are applicable to plates for which the axial load is 
a compression load. More precisely, is given as a 
function of L 2 for selected values of L 3 , and three differ- 
ent curves are shown in each figure that correspond to 
different values of p. The curves in Fig. 28 and 29 are 
for specially orthotropic plates (y = 5 = 0) with values of 
p = 0.5, 1, and 1.5 and with values of p = 2, 2.5, and 3, 
respectively. Curves are given in these Figures that apply 
to both positive and negative shear loads. Similarly, the 
curves in Figs. 30 and 31, and in Figs. 32 and 33, are for 
plates with y = 6 = 0.2 (slight anisotropy) and y = 6 = 0.6 
(substantial anisotropy), respectively, and for P = 2, 2.5, 
and 3. Moreover, the results in Figs. 30 and 32 are pri- 
marily for positive shear loads, and those in Figs. 3 1 and 
33 are for negative shear loads. The distinction between 
positive and negative shear loads, and the differences in 
the corresponding results, is caused by the presence of 
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flexural anisotropy. In each of Figs. 28-33, the curves 
terminate at points on the gray, solid curves given by 

which represent the wide-column buckling 

modes. 

Examples 

To illustrate the use of the generic Figures that have 
been presented herein, First consider a simply supported 
[(±45/90 2 ) g ] s laminate made of IM7/5260 graphite-bis- 
maleimide material. From a laminate analysis code, or a 
set of charts like those shown in Figs. 14-21, one could 

get P = 1.3, y = 5 ~ 0, ft, = 2.27 x 10 6 /°F, a 2 =0.71 x 

10' 6 /°F, and ft 3 = 0. For this laminate, the only destabi- 
lizing loads are biaxial compression loads that are ob- 
tained for uniform heating. For this case, Lj = 1, L 3 = 0, 
and L 2 = 0,27 is obtained from Eq. (47). For L, = 1, K, 

= p CT . The value of is obtained from Fig. 23 by inter- 
polating the results in the Figure for P = 1 and P = 1.5. 
For P = 1 and L 2 = 0.27, Fig. 23 gives K x = 2.9. Similar- 
ly, for P = 1.5 and L 2 = 0.27, Fig. 23 gives K, = 3.5. In- 
terpolating these values gives K* = 3.3 for P = 1.3. 

Substituting K* = 3.3 for p CT and ftj = 2.27 x 10 6 /°F into 

Eq. (49) gives -^jr©o = 1-45 x 10 6o F. Next, let the 

plate width be given by b = 32 in. such that the plate 
width-to-thickness ratio is b/t = 100. For this plate, 

0o « 1 19°F. 

Next, consider a simply supported [±45/90 2 ] s lami- 
nate made of IM7/5260 graphite-bismaleimide material. 
From a laminate analysis code, or a set of charts like 

those shown in Figs, 14-21, one could get P = 2, y= 5 « 

0.2, ft, = 2.39 X 10' 6 /°F, 0-2 = 1.1 1 X 10‘ 6 /°F, and ft 3 = 
0. For this laminate, the only destabilizing loads are also 
biaxial compression loads that are obtained for uniform 
heating. Thus, for this case, L, = 1, L 3 = 0, and L 2 = 0.47 

is obtained from Eq. (47). Again, for = 1, = p CT . 

The value of K x is obtained directly from Fig. 30 and is 
given by = 2.1, which corresponds to a wide-column 

buckling mode. Substituting K x = 2.1 for p CT and ft, = 

2.39 X 10* /°F into Eq. (49) gives ^r 0 o = 0.88 x 10 6 

°F. Next, let the plate width be given such that the plate 
width-to-thickness ratio is b/t = 100. For this plate, 

0" « 72°F. 

Finally, consider a simply supported [(-h45 2 /90 2 ) 8 ] s 
unbalanced laminate made of IM7/5260 graphite-bisma- 


leimide material. From a laminate analysis code, or a set 
of charts like those shown in Figs. 14-21, one could get 

P= 1.3, y= 0.5, 8 = 0.3, ft, = 2.27 x 10 fi /°F, ft 2 =0.71 

x 10 6 /°F, and ft 3 = 0.38 X 10 6 /°F. For this laminate, a 
state of biaxial compression and positive shear loads 
arise for uniform heating. Thus, for this case, L, = 1, L 2 
= 0.27 is obtained from Eq. (47), and L 3 = 0. 17 is ob- 
tained from Eq. (48). Again, for = 1, = p CT . From 

a set of charts like the ones presented herein, one could 

get K x = = 2.6. Substituting = 2.6 for p„ and ft, 

= 2.27 x 10 6 /°F into Eq. (49) gives pjy©o = 1. 15 x 10 6 

°F. Next, let the plate width be given such that the plate 
width-to-thickness ratio is b/t = 100. For this plate, 

©o a; 95°F. Comparing this result with the result for the 
corresponding [(±45/90 2 ) 8 ] s laminate discussed previ- 
ously shows that the unbalanced laminate is less buck- 
ling resistant, as expected because of its higher degree of 
anisotropy and the presence of a shear load. For the case 
of uniform cooling, biaxial tension and negative shear 
loads are induced by the restrained deformation. The 
load factors for this case are obtained from Eqs. (50) and 
(51), where = -1. These equations yield L 2 = -0.27 and 
L 3 = -0. 17. Examination of the load factors indicates that 
the axial and transverse tension loads are approximately 
6 and 1.6 times the magnitude of the destabilizing shear 
load. For this loading condition, the tension loads are 
dominant and elastic buckling is not possible. 

Concluding Remarks 

An analytical approach for synthesizing buckling 
results and behavior for long, balanced and unbalanced 
symmetric laminates that are subjected to uniform heat- 
ing or cooling, and fully restrained thermal expansion 
and contraction has been presented. A nondimensional 
buckling analysis for long flexurally anisotropic plates 
that are subjected to combined loads has been described 
and useful nondimensional parameters have been pre- 
sented. In particular, stiffness-weighted thermal-expan- 
sion parameters have been presented that can be used to 
determine critical temperatures, for a wide range of lam- 
inate constructions, in terms of physically intuitive, well- 
known mechanical buckling coefFicients. Moreover, the 
effects of membrane orthotropy and membrane anisotro- 
py on the mechanical, thermally induced prebuckling 
stress state have been determined. 

A large number of results have been presented 
herein for some common laminates that are intended, to 
some extent, to facilitate a structural designer’s transition 
to the use of the generic buckling design curves that are 
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included in the paper. Many of the results were previ- 
ously unknown. In addition, several results have been 
presented that show the effect of laminate construction 
on the buckling behavior, and several cases are presented 
that indicate when a laminate will buckle because of uni- 
form cooling. Results of this type could be important in 
the design of vehicles that use liquid fuels. Generic 
buckling design curves have also been presented that 
provide physical insight into the buckling problem of the 
present paper in addition to providing useful design data. 
In addition, examples have been presented that demon- 
strate the use of the generic design curves. Overall, the 
analysis approach and generic results that have been pre- 
sented indicate the effects or characteristics of laminate 
thermal expansion, membrane orthotropy and anisotro- 
py, and flexural orthotropy and anisotropy on laminated- 
plate buckling in a very general and unifying manner. 
Although the results are based on classical laminated- 
plate theory and have been demonstrated for infinitely 
long plates, the approach is applicable to more sophisti- 
cated plate theories that incorporate effects such as trans- 
verse-shear flexibility and can be used for Finite-length 
plates. 
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(a) Destabilizing loading system (b) Subcritical loading system 

Fig. 1 Sign convention for positive-valued stress resultants. 



(a) Boundary conditions (b) Positive values for equivalent mechanical loads 


Fig. 2 Mechanical loads in a plate fully restrained against axial thermal expansion or contraction 
caused by uniform heating or cooling. 
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Fig. 3 Stiffness-weighted laminate thermal-expansion parameter Q, for laminates made of IM7/5260 material. 
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Fig. 4 Ratio of stiffness-weighted thermal-expansion parameters and 0 ^ for laminates made of IM7/5260 material. 



Fig. 5 Ratio of sti ffness- weigh ted thermal-expansion parameters and ^ for laminates made of TM7/5260 material. 
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Fig. 6 Stiffness-weighted laminate thermal-expansion parameter &, for balanced and unbalanced laminates with 
angle-plies and unbalanced, unidirectional off-axis laminates made of IM7/5260 material (m = 1, 2, ...). 



Fig. 7 Ratio of stiffness-weighted thermal-expansion parameters fit, and fit, for balanced and unbalanced laminates 
with angle-plies and unbalanced, unidirectional off-axis laminates made of IM7/5260 material (m - 1, 2, ...). 
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Fig. 8 Ratio of stiffness-weighted thermal-expansion parameters ft , and ®3 for [(+0)2 m ] s an d [+02 Z ; O' ; 9O] S unbalanced, 
unidirectional off-axis laminates made of IM7/5260 material ( m = 1, 2, ...). 
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Fig. 9 Critical temperature change for simply supported laminates made of IM7/5260 material, fully 
restrained against thermal expansion and contraction, and subjected to uniform heating. 
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Fig. 10 Critical temperature change for clamped laminates made of IM7/5260 material, fully restrained 
against thermal expansion or contraction, and subjected to uniform heating. 



Fig. 1 1 Critical temperature change for simply supported and clamped angle-ply laminates made of IM7/5260 
material (m = 1, 2, ...). 
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Fig. 12 Critical temperature change for simply supported and clamped laminates with angle-plies made of IM7/5260 
materia] (m = 1, 2, ...). 
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Fig. 13 Critical temperature change versus plate width for simply supported and clamped [(±45/0/90) m ] s 
laminates made of IM7/5260 material, fully restrained against thermal expansion and contraction, 
and subjected to uniform heating. 
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Fig. 14 Effects of lamina material properties on stiffness-weighted laminate thermal-expansion parameters, 
for [(±45/0/90) m ] s laminates. 



Fig. 15 Effects of lamina material properties on parameter ratio ft, / A, for [(±45/0/90) m ] s laminates. 
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Fig. 16 Effects of lamina material properties on stiffness-weighted laminate thermal-expansion parameter S, 
for [(±0) m ] s balanced, angle-ply laminates and [(+0)2 m ] s unbalanced, unidirectional off-axis laminates. 
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Fig. 17 Effects of lamina material properties on stiffness- weighted laminate thermal-expansion parameter 3, 
for [(±0) m l s balanced, angle-ply laminates and [(+0)2 m ] s unbalanced, unidirectional off-axis laminates. 
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Fig. 18 Effects of lamina material properties on stiffness-weighted laminate thermal-expansion parameter Q, 
for [(±0) m ] s balanced, angle-ply laminates and [(+0)2 m ] s unbalanced, unidirectional off-axis laminates. 
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Fig. 19 Effects of lamina material properties on nondimensional flexural anisotropy parameter y for [(+0)2 m ls 
laminates (m = 1,2, ...). 
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Fig. 20 Effects of lamina material properties on nondimensional flexural anisotropy parameter 5 for [(+9)2 m ls 
laminates (m = 1, 2, ...). 
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Fig. 21 Flexural anisotropy parameters for angle-ply laminates made of IM7 / 5260 material. 
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Fig. 22 Effects of orthotropy parameter (3 on buckling interaction curves for specially orthotropic plates (y = 5 = 0) 
subjected to axial compression and transverse tension or compression loads. 
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Fig. 23 Effects of orthotropy parameter (3 and load factor on buckling coefficient for specially orthotropic 
plates (y = 5 = 0) subjected to axial compression and transverse tension or compression loads. 
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Fig. 24 Effects of orthotropy parameter p and load factor on buckling coefficient for specially orthotropic 
plates (y = 5 = 0) subjected to axial compression and transverse tension or compression loads. 
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Fig. 25 Effects of flexural orthotropy parameter P and flexural anisotropy parameters y and 5 on buckling 
coefficients for simply supported plates subjected to axial compression and subcritical transverse tension or 
compression loads. 
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Fig. 26 Effects of flexural orthotropy parameter (3 and flexural anisotropy parameters y and 5 on buckling 
coefficients for simply supported plates subjected to axial compression and subcritical transverse tension or 
compression loads. 


Buckling 
coefficient. * 


^p±2;; 




— Clamped 

— Simply supported 


ft - 

D 12 

+ 2D 6S 

p - 

(Du 

o 2i ) m 

L 2 = 

N„( 



KA 

D22; 

1 = 

N«yl| 

f D u y* 


kA 

iQj. 

tr _ 

(N 

*)J>' 

- 

V(D„ DD 1 ' 4 


Buckling coefficient, K, = x ; ucr ^ 

n'jm; 

Fig. 27 Buckling interaction curves for specially orthotropic plates (y = 5 = 0) with P = 1 subjected to axial 
compression, transverse tension or compression, and shear loads. 


30 

American Institute of Aeronautics and Astronautics 






Buckling 

coefficient, 




Fig. 28 Effects of orthotropy parameter P and load factors on buckling coefficient for simply supported, specially 
orthotropic plates (y = 5 = 0) subjected to axial compression, transverse tension or compression, and shear loads. 
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Fig. 29 Effects of orthotropy parameter P and load factors on buckling coefficient for simply supported, specially 
orthotropic plates (y = 6 = 0) subjected to axial compression, transverse tension or compression, and shear loads. 
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Fig. 30 Effects of orthotropy parameter P and load factors on buckling coefficient for simply supported, specially 
orthotropic plates (y = 5 = 0.2) subjected to axial compression, transverse tension or compression, and shear loads. 



Fig. 31 Effects of orthotropy parameter P and load factors on buckling coefficient for simply supported, specially 
orthotropic plates (y = 5 = 0.2) subjected to axial compression, transverse tension or compression, and shear loads. 
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Fig. 32 Effects of orthotropy parameter (3 and load factors on buckling coefficient for simply supported, specially 
orthotropic plates (y = 5 = 0.6) subjected to axial compression, transverse tension or compression, and shear loads. 
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Fig. 33 Effects of orthotropy parameter [3 and load factors on buckling coefficient for simply supported, specially 
orthotropic plates (7=5 = 0.6) subjected to axial compression, transverse tension or compression, and shear loads. 


33 

American Institute of Aeronautics and Astronautics 




